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Abstract 

We study an interaction of 2D quasiparticles with linear dispersion E = ±u\p\ 
(graphene) with impurity potentials. It is shown that in ID potential well (quantum 
wire) there are discrete levels, corresponding to localized states, whereas in 2D well 
(quantum dot) there are no such states. Scattering cross-section of electrons (holes) 
of graphene by an axially symmetric potential well is found and it is shown that for 
infinetily large energy of incoming particles the cross-section tends to a constant. 
The effective Hamiltonian for a curved quantum wire of graphene is derived and it 
is shown that the corresponding geometric potential cannot form ID bound states. 

1 Introduction 

Monatomic layer of carbon atoms, forming hexagonal lattice (graphene), is studied 
very intensively at present PQ- [3]. "Conical" dispersion law for quasiparticles (the 
name is borrowed from the similar 3D model of a gapless semiconductor) results 
in crucial distinctions of their dynamical characteristics from the corresponding 
characteristics of massive particles. The density of electron states tends linearly to 
zero as a function of energy E, counted from the conical point, that is more rapidly 
than for usual particles in 3D case (v 7 ^). This is a reason to expect that formation 
of bound states in potential wells will be hindered. 
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In the present paper we consider some simple exactly solvable models of ID 
and 2D potential wells from the viewpoint of possibility to form bound states for 
quasiparticles , described by the Hamiltonian 

H = ucrp, (1) 

where er = (o"i,o"2) are Pauli matrices, p = —ihV is momentum operator, u is 
characteristic velocity (for graphene u ~ 10 6 m/sec). It turned out, that quantum 
wire (ID localization) is possible, whereas quantum dot as well as a hydrogen-like 
donor (acceptor) are not possible. ID effective Hamiltonian for a curved wire will 
be derived and it will be shown that geometric potential differs significantly from 
the case of parabolic dispersion law. 

Some problems of interaction of quasiparticles described by the Hamiltonian ([I]) 
with electrostatic potentials have recently been considered. In Ref. [3j the transmis- 
sion coefficient of carriers through ID barrier (p-n junction) is found; same authors 
have shown in Ref. [5] that Friedel oscillations of the charge density around an impu- 
rity atom in graphene differs essentially from the ones in 2D electron systems with 
parabolic dispersion law. Bound states in a symmetric ID potential well have been 
investigated in Ref. [6]. We expound here solution of this problem (together with a 
more general one for an asymmetric well) as a starting point for our consideration 
of a curved quantum wire. 



2 ID potential well 



The motion of electrons in a graphene waveguide, representing 2D stripe with 
straight axis, is described by the equation: 

u(ap)^> + v(y)^> = EV, (2) 

where v(y) is the potential confining the particle in the waveguide (below we suppose 
u = h = 1). Let us look for the solution in the form ^(x, y) = x(y) exp(ipx), where 
x(y) = (Xl>X2) is two-component spinor, whose components satisfy the equations: 



d_ 

dy 



+ P )X2 = (E - v(y))xi 



d_ 

dy 



p)xi = (E- v(y))x2- 



(3) 



After excluding X2i we find: 

d 



dy 



2 +P + 



dy 



ln(E - v(y)) 



0_ 

Oy 



+ P 



Xi = (E - v(y)) xi 



(4) 
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Substituting xi m the form y/~E 
+ p 2 -(E-v(y)) 



v (y)xi> we obtain 



d 2 



dy< 



2 + v"{y)l^ - v'(y)p 



E - v(y) 



Hv'(y)) 2 

A{E-v(y)f 



Xi = (5) 



The equation for function X2 = \J E — v(y)x2 is obtained by means of replacement 
p -p. 

Let us consider the potential v(y) in the form of step function: v (y) = if 
y < —a, v(y) = —vo if \y\ < a, v(y) = —v\ if y > a. In each region the equation © 
is reduced to 

-d 2 X i/dy 2 = [(E - Vi f - p 2 ] X i, 

where Vi = if y < —a, V{ = — i>o,i if |y| < a and y > a correspondingly. Solutions 
of the last equation decreasing when y — » ±oo, have the form xi = Aie^ if y < —a, 
Xi = Bi s'mqy + E>2 cos qy if \y\ < a, xi = A^e^^ if y > a. Here x = \J p 1 — E 2 , 
q = y/(E + v ) 2 -p 2 , Xi = y/p 2 - {E + Vl) 2 . 

Matching conditions Xi\y=±a-o = Xi\y=±a+o, i = 1, 2 lead to the following equa- 
tion, determining the spectrum E = E u (p), where v is the number of subband of 
transversal quantization: 
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K + p Xl — p 



+ 



P 



Xi — p x + p p 



E E + vi E + v E + vi 



Xi 



E E + vq 



sin(2ga) 



E + v 



p X + p 

+ 



cos(2ga) = 0. 



(6) 



E + vi E 

In the case of symmetric well v\ = 0, xi = x; equation ([6]) is simplified: 

[E(E + vo) - p 2 ] sin(2ga) - xg cos(2ga) = 0. (7) 

Branches E u (p) for symmetric well are shown on Fig.l. Let us note that equation 
([7]) always has the solution E(jp) = \p\ — vq, i.e. q = 0. However this branch of the 
spectrum is not physical since it corresponds to the wavefunction identically equaled 
to zero. 

Investigation of the equation ([7]) leads to the condition of reality of q, which, 
together with obvious condition of reality of x, determines the region, occupied by 
discrete spectrum max{— \p\, \p\ — vq} < E < \p\. From this, in particular, it follows 
that E > —vq/2. It is easy to calculate energy in the following limiting cases: 



1) \p\ < v ,cos(2t> a) / : E (p) 

2) \p\ < v ,cos(2v a) = : E(p) 

3) wo^O a^O: E = \ 



|p|sign{sin(2voa)} cos(2woa), (8) 
P 2 /v , (9) 
{l-2v 2 a 2 + ...). (10) 



In the shallow asymmetric well there are no bound states at any p. There are 
no bound states also in an asymmetric well of any depth at small enough p. 
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Figure 1: Bound states E(p) in the ID symmetric rectangular well; v^a = 5. Dotted line 
shows the bounds of the area, occupied by discrete spectrum. 

3 Effective ID equation in adiabatically 
curved stripe 

Let us consider a stripe with curved axis defined by the equation r = R(x); we 
suppose the stripe to be curved in its plane. Here R(x) = (Ri(x), i?2(x)) is a 
smooth vector-function, x is the natural parameter on the axis (length counted out 
from some fixed point), \dR/dx\ = 1. Let n(x) = (ni(x), ri2(x)) be the unitary 
vector normal to dR/dx. In the neighborhood of axis the curvilinear coordinates 
x, y, defined by the equality r = R(x) + yn(x) can be introduced. The equation 
(El) in curvilinear coordinates has the form: 



(crR') . d 1 . d ik(an) 

= (-t— )-==== + (<rn)(-t— ) - n ./ /. + v{y) 



y/\ — ky dx y/1 — ky dy 2(1 — ky) 



ty = E$>, (11) 



where k is the curvature of the axis of waveguide at the point x, \P = (1 — ky) 1 / 2 ^. 
Let us take into account now that the stripe is curved adiabatically, i.e. formally 
k — > 0. Following the idea [7], we look for the solution of the last equation in 
the form ^ v {x,y) = ^^(x), x u = (X11X2) 1S two-component vector-operator and 
scalar wavefunction ijj u {x) satisfies effective Schrodinger equation 

L u V(x) = Ei> u (x). (12) 
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Now we expand the operator L in powers of curvature L v = + L\ + . . .. It will 
be shown that operator H^ s = H^ s (—id/dx) does not depend on x. It follows from 
[8] that Hg S (p) is an eigenvalue of the following problem: 



d 

{aR')p+{an)(-i—) + v{y) 



where p is a parameter (c-number) 
(712 — in±)x2 we reduce (fT3|) to Q: 



X" 



x"(p, x, y) = H" s (p)x u (p, x, y), (13) 
= (XiiX^)- Using the replacement X \ = 



d 
dy 



X-2 



d 
dy 



Xi = {Hesip) ~ v{y))x2- (14) 



It follows from the last equation that H" s (p) = E u (p), i.e. H^ s does not depend on 
x. Let us choose functions \\, X2 to be real. Then from general formulae [8] and 
the relation {{x v ) + (crn) X v ) y = it follows that 

L<{=[k(x),{2K(p)p + Q(p))} + , K(p) = (xZy%) y , Q(p) = (XiX2> tf ■ (15) 

Here p = —id/dx, [-, •]+ means anticommutator, and (■) means integration over y. 
For even potential K{p) = 0. 

Effective equation for rectangular well at small \p\. For symmetric 
rectangular well at small \p\ one can use dispersion relation ([8]). Then we find 
Q(0) = sign{sin(2uoa)}sign{p} cos(2t;oa)/2. Effective ID equation for small mo- 
mentums takes the form: 

/ Q k(x)\ 

sign{sin(2t;oa)}sign{j5} I — i— H — tp = Etp (16) 

\ ox 2 J 

Thus, the geometric potential for quasiparticles has the form k(x)/2 and, how one 
can see from (I16p . it cannotform bound states. 



4 Axially symmetric potential 

Let us consider the equation ([2]) for the axially-symmetric potential. In the cylindric 
coordinates x = rcoscp, y = rsintp it has the form 

e ~ lV (-''TT " ~ir) = u(r))*i, e* (-£ + ~) * 1 = (E- u(r))* 2 .(17) 
\ Or r oip J \ or r dip J 

We find the solution in the form ^1 = e inip xii ^2 = e^ n+1 ^X2 ; where Xl> X2 satisfy 
equations 

/ d i(n + l)\ ,„ , ,. / d in\ ... . , 

[~ l d~r " ^~^J X2 = (E- v(r)) Xl , + -j Xl = ( E - v(r)) X2 . (18) 
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4.1 Rectangular well 

Let us consider axially symmetric well with constant width: v(r) = —vq if r < a 
and v(r) = if r > a. Excluding \2 from (|18|) . we obtain: 

d 2 1 d n 2 1 

^--7T + — -(^ + «o) 2 ^Xi = 0, r<a, (19) 

These equations have solutions J n (\E + vo\r), N n (\E + vo\r) and J n (\E\r), N n (\E\r) 
correspondingly. The solution regular in the point r = has the form Xi( r < a ) = 
C n J n ((E + v )r). 

The equation (|20p contains only the square of the energy. Thus its solutions 
don't depend on the sign of E. These solutions are equivalent to scattering states 
of the usual radial Schrddinger equation with E > 0. From this circumstance one 
can conclude that there are no bound states in such a potential well. We would like 
to emphasize that this conclusion does not depend on the depth and width of the 
well, i.e. 2D localization of quasiparticals in graphene (quantum dot) is principally 
impossible (naturally, this statement only relates to the region of momenta where 
the Hamiltonian ([2]) and linear dispersion law are valid). The same is true, obviously, 
for any potential decreasing at the infinity, hence, hydrogen-like states of donors or 
acceptors in graphene do not exist. 

Let us consider now the scattering problem. Let the wave with positive energy 
comes from the infinity along the x-axis. Then the wavefunction \I/ far from the 
origin has the form: 

Using the expansion of the plane wave exp(i\E\r cos ip) = 

(1/2) Y^=-oo i n \Hn {\E\r) + Hn\\E\r)\ exp(my), we can represent the solu- 
tion in the form 

* = I V (i n e in nsk 1) Hi l \\E\r) + 2-V 2 Hk 2 \^ l , „ , 
2 n f-U™e^(5i 2) ^ 1) (| J E|r) + 2- 1 /2 jff (2) r - • ' 




Hence 



h { H> ) = ±— £ e^(SW - 2-V 2 ), 
sj 2ir\E\ 



■ — OO 

oo 



-in/4 

M<p) = -_ e ^(^ 2) - 2_1/2 ) ( 23 ) 

V 27T \ h \ n=-oo 
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(2) 

Withfunctional relations for Hankel functions, it is easy to show that S'^_j_ 1 = 

Sn \ Hence ^(v?) = e *^/i(y)- Continuity conditions for the wavefunction at r = a 
lead to the equations: 

S^H^\\E\a)+2~ l l 2 H {2 \\E\a) = C n J n (\E + v \a), (24) 
S^H^QEla) + 2- 1 / 2 H^ l (\E\a) = C n J n+1 (\E + v \a) (25) 

From these equations one can find Sn \ Sn \ Wavefunction of the scattered particles 
has the form ^ ou t — e i ^ r r~ 1 ' 2 (fi ((p), f2(p)), and the current density is given by 
3 out = ^tut^^out = 2|/i((/?)| 2 r -1 {cos ip, sin ip} (we write Cartesian components of 
the current in braces). Thus, the differential cross-section is da /dip = 2\fi(ip)\ 2 . 
Total cross-section is 

2 



8 oo 

\E\ ^ 

1 1 n=0 



<r = 2 / |/ 1 ( ¥J )p^= £ + l 

^ ' ' n=-oo 

J n (|-E + v \a)J n+1 (\E\a) - J n +i(\E + v \a)J n (\E\a) 



J n (\E + v \a)H%{\E\a)- J n+l {\E + vv\a)H { n\\E\a) 



(26) 



Similar calculations lead to the same formula for the total cross-section for particles 
with negative energy. In the low-energy limit \E\a <C 1 the asymptotic formulae 
J (z) = l + 0(z 2 ), Ji{z) = z/2 + 0(z 2 ), N (z) = const + 21n(z)/7r + 0(z), Ni(z) = 
-2/(irz) +O(z ) result in 

a ^ 2 (^4) 2 \E\a 2 (27) 
V Jo{v a) ) 

In Fig. 2 we show the scattering cross-section as a function of \E\a. One can see 
that cross-section has resonances. We will study below the opposite limiting case 
I -El a 3> 1 for which a different method is more convenient. 



4.2 Green function and integral equation for the scat- 
tering problem 

Green function of the operator ([I]) is 2 x 2-matrix G(r). It satisfies the equation 
[(ap)-E]G(r-E)=5(r)I, G = ^) , I = (j J) . (28) 

Solving this equation we find 

rWnfPi^ I w.ic—if ( wW V 



G(r . E) = I ( iEHq (\E\r) \E\e-^{H^)\\E\r) 
4 \\E\e^ {H^ ] )' {\E\r) %EH§\\E\r) 



(29) 




Figure 2: Scattering cross-section by rectangular axially symmetric well with radius a, 
v a = ir. 



8 



For large \E\r Green function has the asymptotic form 



g( ^)= v h, (I ^»(^) ^j E) y ( 3o, 

The scattering is described by the integral equation 

#(r; E) = * in (r; E) - J G(r - r'; E)v{r')^{r'; E)d 2 r' , (31) 

where *$>i n is the wavefunction of incoming particles. In the first Born approximation 
we obtain: 



v a jir\E\ jmr+n/i) ( 1 + e 



*(r;E) = VUr-.E) - ^-L^WV ^ - + ^ J * („), (32) 

where g = 2|i£| sin(y>/2). Scattering cross-section is 

der 27r\E\v 2 ,a 2 9/ t9/ , nv^a 2 9 , , . <,,„, „. , , , N 
■ cos 2 (<^/2) Jf(qa) = —^—cot 2 (<p/2)J 2 (2\E\asin((p/2)), 



dip q 2 yr ' ' 1Vf ' 2\E\ 

o = f^B. ^[i/a, 3/2; 2, 2, 3; -4|£|V], (33) 



where 2F3 is the generalized hypergeometric function. 

When \E\a — ► 00 scattering cross-section tends to the constant cr — ► 16t; 2 a 3 /3. 
Note that for usual particles the Born scattering cross-section by a short-range 
potential tends to zero when energy increases [9], however if one formally supposes 
the mass to be proportional to momentum (to obtain the linear dispersion law), then 
limiting value of the Born cross-section at \E\ — * 00 also turns out to be constant. 
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